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ABSTRACT 

We consider the asymptotic behaviour of the Gauss hypergeometric function when several 
of the parameters a, b, c are large. We indicate which cases are of interest for orthogonal 
polynomials (Jacobi, but also Krawtchouk, Meixner, etc.), which results are already available 
and which cases need more attention. We also consider a few examples of 3F2 functions of 
unit argument, to explain which difficulties arise in these cases, when standard integrals or 
differential equations are not available. 
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1. Introduction 

The Gauss hypergeometric function 

2 F 1 f a '^^=l + ^+ a(a + l) y+ 1) ^ 2 +•^^=f:W ^ %^ 

V c J c c(c + 1)2! ^ (c)„ n\ 

where 

. , T(a + n) , .„ , (— a\ 

<°>" = J TM i = ( - 1 > "'(»)' <L2) 

is defined for \z\ < 1 and c / 0, —1, —2, .... 
A simplification occurs when b = c: 

a F X r°^*) =(!-*)--• (i-3) 

When a or b are non-positive integers the series in (1.1) will terminate, and F reduces to 
a polynomial. We have 

F f-n,b ) = y (z^k% ^ = y { - ir ( n )^n zm (L4) 

V c J (c) m m\ ^ \mj(c) m 
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The value at z = 1 is defined when 5ft (c — a — 6) > and is given by 

/o,b \_ r(c)r(c-a-6) 

2 1 U ' 7 " r(c - a)r( c - 6) • (L5) 

For the polynomial case we have 

,Fj- n ; b ;l)=^^, n = 0,l > 2 > ... . (1.6) 

Generalizing, let p, g = 0, 1, 2, . . . with p < q + 1. Then 

fai,...,a p _ \ _ (ai)n • • • (« P )n ^" n 

This series converges for all z if p < g + 1 and for \z\ < 1 if p = q + 1. 

Sometimes we know the value of a terminating function at z = 1 (Saalschiitz's theorem): 



3 r 2 



-i, a, 6 . \ = (c- a) n (c- b) n , . 

a, l + a + 6-c-n ' ^ ( c ) n ( c - a - 6) n ' 1 j 



where n = 0, 1, 2, . . . . 

The behaviour of the Gauss hyper geometric function F(a, b; c; z) for large \z\ follows from 
the transformation formula 



2-f 1 



a,b_ \ _ r(c)r(6 - a) /a,a-c+l.l 
c ' ^ " r(6)r(c-o) 1 J 2 ^Va-6+l'z 



(1.9) 



(1.10) 



, T(c)T(a-b) , b /6,6-c+l. 1 
+ T(a)r( C -6) 2Fl V 6-a + l ' * 

where |ph(— z)| < n, and from the expansion 

\ c y c c(c + 1) 2! 

with 

(a)n = E r^r = ( " irn! (r)' n =°. i .2.-. ( Ln ) 

and where |z| < 1. We always assume that c / 0, —1, —2, .... 

The asymptotic behaviour of F{a, b; c; z) for the case that one or more of the parameters 
a, b or c are large is more complicated (except when only c is large). Several contributions 
in the literature are available for the asymptotic expansions of functions of the type 

F(a + eiA, b + e 2 A; c + e 3 A; z), ej=0,±l, A — > oo. (1-12) 
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In this paper we give an overview of these results, and we show how this set of 26 cases 
can be reduced by using several types of transformation formulas. 

In particular we indicate which cases are of interest for orthogonal polynomials that can be 
expressed in terms of the Gauss hypergeometric function. The Jacobi polynomial P^ 01 '^ (x) 
is an important example, and we mention several cases of this polynomial in which n is a 
large integer and a and/or (3 are large. When a or [3 are negative the zeros of Pi a '^(x) may 
be outside the interval [—1, 1], and for several cases we give the distribution of the zeros in 
the complex plane. 

We mention a few recent papers where certain uniform expansions of the Gauss func- 
tion are given (in terms of Bessel functions and Airy functions), from which expansions the 
distribution of the zeros can be obtained. Many cases need further investigations. 

Of a completely different nature is the asymptotics of generalised hypergeometric func- 
tions. We consider a few cases of 3-F2 terminating functions of argument 1 and —1, and show 
for these cases some ad hoc methods. In general, no standard methods based on integrals or 
differential equations exist for these quantities. 



2. Asymptotics: a first example 

Consider the asymptotics of 



, a, j3 + A , 
2-Fi I 7 + A 5 z ) , A -> oo. 



(2.1) 



We use (3 + A ~ 7 + A and we try (using (1.3)) 



„ . a,(3 + X \ „ f a,P + X 1 . „ 9 

2F 1 [ 7 + A ,zl~ 2 Fil p + x ,z)-(l-z) ? 



(2.2) 



Observe first that, if 3? (7 — a — (3) > 0, 

a, (3 + A_ 



, , 1 i = r(7 + A)r( 7 -a-/?) 

2 1 7 + A ' J r( 7 + A - o)r( 7 - f3) ' 



(2.3) 



We see that (2.1) cannot hold for z ~ 1 (if ^(7 -a-/3)>0). 
However, we can use 



2 F l{ a l b ;z)=(l-z)-« 2 Fj a > C - b - Z 



z-l 



(2.4) 



This gives 



. a, 13 + X \ _ , / a, 7 - , z 

2Fl{ 7 + A ^M 1 "*) 2i?1 ( 7 + A 



(l-z)- a 
(l-z)- a 



1 + a( 7 -/?) z + (o)(a + 1)( 7 - /3)( 7 - /3 + 1) 



7+A z-l 



(7 + A)( 7 + A + l)2! 



z - 1 



(2.5) 



l + ^^^- + 0(A- 2 ) 
7 + A z — 1 
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as A — > oo with z fixed. This is the beginning of a complete asymptotic expansion, which 
converges if \z/(z — 1)| < 1. It is an asymptotic expansion for A large, and all fixed z, z / 1. 
We can also use other transformation formulae: 

2Fl ( a * ; z) = (1 - z)-\ Fl ( C 7' 6 ; -± T ) = (1 - z) c ~ a ~ b 2 F 1 ( C " ^ ~ & 5 *) 
These give, with (2.4), the three relations 

fa,(3 + X \_ n x-q F / 'q,7-£. z \ 

= (1 - ,)-^ ( 7 + A " + Y + A ; j^j) = (1 - ( 7 + A - P ; ,) 

(2.6) 

We see in (2.7) that the large parameter A can be distributed over other parameter places. 
In the present case only the second form is suitable for giving an asymptotic expansion when 
using the Gauss series. The third form gives a useless Gauss series (for large A). In the final 
form the Gauss series converges at z = 1 if 9ft (a + (5 — 7) > 0, but the series does not have an 
asymptotic property for large values of A. 

The transformation formulae in (2.7) are an important tool for obtaining asymptotic 
expansions. We will use these and other formulae for investigating all cases of (1.12). 

So far, we mention the Gauss series for obtaining an asymptotic expansion. However, to 
obtain an asymptotic expansion of (2.1) that holds uniformly if z is close to 1, we need a 
different approach, as will be explained in later sections. 

3. Some history and recent activities 

[1] Watson (1918) studied the cases 

/ a + X ,p + \ 2\ ( a + X ,[3-\_l-z\ ( a ,b \ 

All results are of Poincare type (i.e., negative powers of the large parameter) and hold 
for large domains of complex parameters and argument. Watson began with contour 
integrals and evaluated them by the method of steepest descent. 
[2] Luke (1969, Vol. I, p. 235) summarizes these results, and gives many other results, 
also for higher p F q — functions. In particular Luke investigates "extended Jacobi poly- 
nomials" , which are of the type 

f -n,n + X,ai,...,a p \ 

p+2i H &!>■■■ A 5 7' 

For p = 0, q = 1 and integer n this is the Jacobi polynomial. 
[3] Jones (2001) considers a uniform expansion of 

/ a + X ,(3-\l 1\ 
2Fl { c ' 2 " 2 Z ) 
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and gives a complete asymptotic expansion in term of I— Bessel functions with error 
bounds. Applications are discussed for Legendre functions. 
[4] Olde Daalhuis (2001) and (2002) give new uniform expansions of 

F ( a ^~ X . 7 \ 
2i H 7 + A '~ z ) 

in terms of parabolic cylinder functions, and of 

,*( + *f + ";-.) 

in terms of Airy functions. 

[5] Other recent contributions on uniform expansions for Gauss functions follow from 
Ursell (1984) (Legendre functions), Boyd & Dunster (1986) (Legendre functions), 
Temme (1986) (see later section), Frenzen (1990) (Legendre functions), Dunster 
(1991) (conical functions), Dunster (1999) (Jacobi and Gegenbauer polynomials), 
Wong and co-workers (1992) (Jacobi polynomials), (1996) (Pollaczek), (1997) (Ja- 
cobi function), (1998) and (1999) (Meixner), (2000) (Krawtchouk) , (2001) (Meixner- 
Pollaczek) . 

[6] On Bessel polynomials, which are of 2-Fo - type, see Wong & Zhang (1997) and 
Dunster (2001). On Charlier polynomials ( 2 F -type) see Bo Rui & Wong (1994). 
On Laguerre polynomials (i-Fi— type) see Frenzen & Wong (1988) and Temme 
(1990). 
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Table 1. All 27 cases of 2 Fi(a + eiA, b + e 2 A; c + e 3 A; z), ej = 0, ±1 . In the 
table we write ej = 0, ± . 



4. From 27 cases to only a few 

Skipping the dummy ( ), using the symmetry between the a and b parameters: 

2 F 1 ( a l b ; Z )= 2 F 1 (\ a ; Z ), (4.1) 



c I \ c 
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and by using the transformation formulae of (2.6) we reduce the 27 cases to the 8 cases: 
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This set is obtained by using (4.1) and (2.6), that is, only 4 of Kummer's 24 solutions of 
the hypergeometric differential equation. 

A further reduction of the set of 8 cases can be obtained by other solutions of the differ- 
ential equation. 

Using the connection formula 

f a ,b \_ r( c -i)r(q- c + i)r(b- c +i) ! b /i- a ,i-6 ' 

2Fl {c> z )- r(a)r(6)r(i - c) z {l ~ z) 2Fl { 2-c > z 

T(b-c+l)T(a-c+l) f a,b 
r(a + 6-c+l)r(l-c) 2 1 Va + 6-c+l' 

we see that the second case (0 — ) in fact reduces to ( + ). 
Similarly, we have 

jp fa,b \ j T(c)T(b-c+l) { a }a - c+ i 1 

2Fl Vc'7" T(a + b-c+inc-a) Z \a + b-c+V z 



(4.2) 



(4.3) 



r(o)r(6-o + l) v ; V b-a+l 

which reduces the third case (0 + ) again to ( + ). 
Finally, we use 

F \_ rC ^ nc-mb-c+l)T(l-a) , b /l-o,l-6. \ 

2jFi U'J" r(6)r( c -a)r(i- c ) — z {1 ~ z) 2Fl { 2-c > z ) 

, Jc-a)ni r ( l - a ) T (b-C+l) ^c-a-b p / 1 - a, C - _ l\ 

+ C r(l-c)r(6-a + l) Z 2Fl V 6-a + l 'ij' 

(4.4) 

which reduces the forth case (0 + — ) to the fifth case (0 — + ). 
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Consequently we have the five remaining cases: 
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These five cannot be reduced further by using relations between the functions 

/ a + e 1 X,b + e 2 A , 
^ I. c + e 3 A ;Zi '-.-»-±l- 

However (Olde Daalhuis (2001)), using 



2 F l( a ' 6 ; , 



r a^ r (c) r ( fc -«) -a p /a,a-c + l. 1 
r(6)r(c-o) 2 1 V a- 6 + 1 ' z 



bm T(c)T(a-b) _ b 
+ e T(a)T(c-b) Z 2Fl 



6,6-c+l 1 
b-a+1 ' z. 1 ' 



and 



f a, & . ^ _ Jc-b) ^(c)r(b - a) b _ _ b ( 1 - &, c - b . 1 

21 V C / r(6)r(c-a) 1 J ^a-fe+l'z 



T(a)T(c - b) 



z a ~ c {\ - z) c - a - b 2 F 1 



1 — a,c — a 1 
b-a+1 'I 



(4.5) 



(4.6) 



we see that the cases ( + + — ) and ( — — + ) both can be handled by the case 
(+ 2+ 0). 

Hence, the smallest set is the following set of four: 
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if we accept the special case D: 



,. . a + X,b + 2X 
2*1 ( c \z 



8 



4.1. Quadratic transformations 

For special combinations of the parameters a,b,c other types of relations between 2-fi- 
functions exist. For example, we have the quadratic transformation 

O.6. \ _ n .-a/2 „ / i a ' 6 - i°. ^ A 



26 '7 v ; z 1 ^ b+\ ' Az-AJ" 
For large a the left-hand side is of type 

( + ) = ( + ) = ( + ) = A, 
whereas the right-hand side is of type C = ( + — ). 
5. Which special functions are involved? 

The next step is to indicate which special functions are associated with the above four cases. 

5.1. Legendre functions 

We have the Legendre function 

Case A: (00 + ) 

We see that case A occurs when \i — > —00. Several relations between the Legendre 
function and the Gauss function give also the case fi — > +00. 

Case B: ( - + ) 

Because 

(- + -) = (0 + -) = (0 - +) = B 

(see (3.5)), we see that case B occurs when in P^(z) both parameters v and \x tend to +00. 
Case C: ( + - ) = ( - + ) 
This also occurs for P£(z), with v large and \i fixed. 

In the results for Legendre functions in the literature rather flexible conditions on the 
parameters jj, and v are allowed. For example, in Dunster (1991) it is assumed that the 
ratio v I [i is bounded. 

Any hypergeometric function, for which a quadratic transformation exists, can be ex- 
pressed in terms of Legendre functions. So, many special cases, and mixed versions of the 
cases A,B,C,D are possible for Legendre functions. 

5.2. Hypergeometric polynomials 

Several orthogonal polynomials of hypergeometric type (see the Askey scheme in Koekoek 
& SWARTOUW (1998)) have representations with Gauss hypergeometric functions. Next to 
the special cases of the Jacobi polynomials (Gegenbauer, Legendre) we have the following 
cases. 
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Pollaczek 



P n (cos 9; a, b) = e^F, ( " n ' 3 + ^ ; 1 - e" 2 ^) 



.... a cos 9 + b 
(f)(9) = 



sm( 



For uniform expansions in terms of Airy functions see Bo Rui & Wong (1996). The asymp- 
totics is for n — ► oo, = t/y/n with t bounded, and bounded away from zero. 



Meixner-Pollaczek 



For uniform expansions in terms of parabolic cylinder functions see Li & Wong (2001). The 
asymptotics is for n — ► oo, x = an with a in a compact interval containing the point a = 0. 



Meixner 



M n (x;(3,c)= 2 F 1 { n ' X ;l-- 



c 



In Jin & Wong (1998) two uniform expansions in terms of parabolic cylinder functions are 
given. The asymptotics is for n — > oo,x = an with a > in compact intervals containing 
the point a = 1. 



Krawtchouk 



K n (x;p,N) = 2 Fi ( "' iv a; ;^)> n = 0,l,2,...,JV. 

A uniform expansion in terms of parabolic cylinder functions is given in Li & Wong (2000). 
The asymptotics is for n — ► oo, x = AiV with A and ^ = n/iV in compact intervals of (0, 1). 

5.3. Jacobi polynomials 

The Jacobi polynomial has the representation: 

*<"»<*) - (" °:^ B+i i |d-*)). (") 

and we see that case C: (+ — 0) = (— + 0) applies if n — > oo. Uniform expansions 
use Bessel functions for describing the asymptotic behaviour at the points x = ±1 (Hilb-type 
formula, see Szego (1975) or for complete asymptotic expansions Wong & Zhang (1996)). 
A well-known limit is 



lim a -n/2p(a+a,b+a) 



_ H n (x) 
y/a~J 2 n n! ' 
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where H n (x) is the Hermite polynomial, a special case of the parabolic cylinder function. 
Approximations of H n (x) are available in terms of Airy functions. 

Several other limits are known (see Koekoek & Swartouw (1998)). In Chen & Ismail 
(1991) and Dunster (1999) asymptotic expansions are given for large positive parameters 
a and/or (3. 

Of particular interest are asymptotic expansions for large n and large negative values of 
the parameters a, (3. See Bosbach L Gawronski (1999) and Marti'nez-Finkelshtein 

ET AL. (2000). 

With a = (3 = —n we have in (5.2) the case 

( ) = ( -) = ( + ) = A 

These are non-classical values for the parameters. We have, taking a = (3 = — \ — n, 

This representation shows that all zeros are on the imaginary axis. 
Other representations can be used: 

^>-(T)(^)~W a+n «V~"^H 

^(":')(^)""(-w(' + ^r- ; ^4 

and we see that case D : ( + 2 + ) applies if n — ► oo and, in the first case, a is 
and (3 equals b — 3n, and similarly in the second case. 

See Figure 1 for the zero distribution of the Jacobi polynomial 

Pi a '^{z), n = 30, a = \, /3 = -3n+l. 



Ln/2j 

(-ir E 



n\n\1~ n ~ 2 m T- n— 2m 



m=0 



m\m\(n — 2m)! 



i » t t • • • • 



1 5 10 



Figure 1. Zeros of P$? ,p) {z), n = 30,a = ^, /3 = -3n + l. 
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The case B : ( — +) = ( — 0+) applies in the representations 
P^\x) = 



a + (3 + 2n\ fx -I 

n J 
a + f3 + 2n\ (x + l 
n 



2 Pi 
2 Pi 



-n, —a — n 



-2n- a- /?' \ -x 
-n, —(5 — n 2 



-2n-a-f3' 1 + . 



(5.4) 



when a = —a — n,(3 = — b — 2n (first case) or @ = —b — n,a = —a — 2n (second case), 
with, again, non-classical parameter values. Olde Daalhuis (2002) gives an expansion of 
the Gauss function for this case in terms of parabolic cylinder functions. 
See Figure 2 for the zero distribution of the Jacobi polynomial 



P<?'®(z), n = 25, 



a 



-n+i, f3 = -2n + l. 



As n — > 00, the zeros approach the curve defined by 



1. 



0.5 



-0.5 




Figure 2. Zeros of P^'^z), n = 25,a = -n+^, /3 = -2n + l. 



Uniform expansions for Jacobi polynomials are available in terms of J— Bessel functions. 
Details on uniform expansions of Jacobi polynomials in terms of Bessel functions can be 
found in Wong & Zhang (1996), Dunster (1999), and in cited references in these papers. 

5.4. Gegenbauer polynomials 

For Gegenbauer (ultraspherical) polynomials we have several representations. For example, 

r A^_ (2A)n f-n,n + 2X l- x \ 

and for A = — n, n — > 00, we have 

( ) = ( -) = ( + ) = A 

We also have 

^w = (-ir^(-"'r A ^ 2 )' 
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and for A = — 2n, n — > oo, we have 

( _ _ o ) = ( + - o ) = a 

Also in this case a "sort of" quadratic transformation is used. 
5.5. Asymptotics of the A— case 

We give a few steps in deriving the expansion for the case A = ( + ). As always, we can 
use the differential equation or an integral representation. We take 

„ (a,b \ Tjc) /■' «>-'(!- if-*- 1 

where JRc > !R& > 0, |ph(l — z)\ < ir. A simple transformation gives the Laplace transform 
representation 

with A — > oo and where 

6-1 



/(*) = ( f T i ) " e (1 - c)t (i-^ + ^)" a - 



This standard form in asymptotics can be expanded by using Watson's lemma (see Olver 
(1974) or Wong (1989)). For fixed a, b, c and z: 

a,b _\ T(c + A) ^ f (p) s 



2i?1 { c + A ; Z ) ~ T(c + A - 6) ^ /s( " ; A 6+s ' 
where /o(^) = 1 and other coefficients follow from 

oo 

/(*) = ^2fs(z)t 8 , \t\ <min(27r,|io|), i = m^-. 



s=0 



When z is large, the singularity t is close to the origin, and the expansion becomes 
useless. 

We can write 

F ( r ( c+A ) rt^M* 

^{c + X^ Z )-T(b)T(c + X-b)J (t + Cr ' 

where 

C = -io = ln^, s(t) = (t + C)°/(t), 
with (7 regular at t = t . Expanding g(t) = 9s(z)t s gives 

^ ( Xx ! z ) ~ r(c + A-6) S 5s(^) (fo) s C 6 " a+ ^(& + 5 - a + 1 + S , CA), 
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in which U is the confluent hypergeometric function. This expansion holds if A — > oo, 
uniformly with respect to small values of £, that is, large values of z. 

This method can be used to obtain an expansion of the Legendre function P^ m (z) for 
large values of m; see Gil, Segura & Temme (2000), where we used this type of expansion 
for computing this Legendre function for large m and z. 

Other sources with related expansions for this type of Legendre functions are Boyd & 
Dunster (1986), Dunstger (1991), Frenzen (1990), Olver (1974) and Ursell (1984). 

5.6. Another A— case 

The case ( — ) is very important for all kinds of orthogonal polynomials and can be 
reduced to the previous case A = ( + ) . It is of interest to give a direct approach of 



2^1 



-n, b _ 



n — > oo, 



(5.6) 



whether or not n is an integer. First we recall the limit 

a, b 



lim 2-^1 

b— >oo 



; z/b] = iFx 



; z 



which may be used as a definition for the \F\— function (the Kummer or confluent hyperge- 
ometric function). We are interested in the asymptotics behind this limit, and we expect a 
role of the Kummer function when in (5.6) n becomes large, and z is small. We again can 
use (5.5), and consider 

I n = f 1*- l (l-t) c - b - l (l-tz) n dt, 
Jo 

where Kc > > 0, |ph(l — z)\ < it. We transform 

\-tz = {1-zY 



= uln(l-z) 



Then we have 



(1-z) 



c-b-l 



ln(l - z) 



—z 



c-1 



J n = f f{u)u b -\l-u) 
Jo 



/(«) 



1 - (1 - z) u 
— nln(l — z) 



n 6-1 r 



i-a 



\u-l 



-i c-b-l 



(1 -u) ln(l - z)_ 



u) = (n+ l)ln(l - z). 
The function f(u) is holomorphic in a neighborhood of [0, 1]. Singularities occur at 

2kiri 2miri 



In(l-z)' 



1 + 



ln(l - z) 



, fc,mG2Z\{l}. 



So, when z ranges through compacta of (D \ {1}, the singularities of / are bounded away from 
[0,1]. If $tuj < the dominant point in the integral is u = 0; if ^tuo > 0, then u = 1 is the 
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dominant point. To obtain an asymptotic expansion for large n that combines both cases, 
and which will give a uniform expansion for all z, \z — 1| > <5 > 0, / should be expanded at 
both end-points and 1. 

More details can be found in Temme (1986), where the results have been applied to a 
class of polynomials biorthogonal on the unit circle. 

5.7. Asymptotics of the B— case 

We consider a simple case: 



( n + 2 ! -*) = ( n + X ) jf [(!-*)(! + *)]' 



dt, 



with z near the point 1. 
We write 



[{i-t){i + zt)]' 



-n<f>(t) 



and we have 



b'(t) 



, cj>(t) = -\n[(l-t)(l + zt)\, 



-2tz -1 + z 



(5.7) 



{i-t)(i + zty 

We see that the integrand has a peak value at to = 

So, if z = 1 the peak is at t = 0, if z > 1 then to G (0, 1), and if z < 1 then t < 0. See 
Figure 3. 




z = l/4 




z=l 




Figure 3. (1 — t)(l + zt) for a few values of z. 



The same situation occurs for the integral 

e - n{ i w *- aw) dw 



f 

Jo 



(5.8) 



which has a peak value at wo = a. This integral is an error function (a special case of the 
parabolic cylinder functions). 

In a uniform expansion of (5.7) that holds for large n and z in a neighborhood of 1 we 
need an error function. This explains why in Case B (in a more general form than (5.7)) 
parabolic cylinder functions may occur for special values of z. 
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We transform the integral in (5.7) by writing 



(1 - t)(l + zt) = e~* 



w +otw 



with the conditions 



t = w = 0, t = 1 if = oo, t = to u> = a. 
The quantity a follows from satisfying the matching of to with a: 

(l-t )(l + zt ) =e-^ w ° +aw °, 



giving 



We obtain 



±a 2 = -ln 



z- 1 

ITT 



sign(a) = sign(z — 1). 



n + 2' 



(n + l) / 



-n(^w 2 — aw) 



f(w) dw, 



where 



(it _ w - a (1 - t)(l + zt) 



diu t — to 2z 
A first approximation follows by replacing f(w) by 

l + z 



/(«) 



2^' 



giving 



(5.9) 



2Fl (n + 2' _Z ) ^ nf ^ J e~ n ^ w2 ~ aw) dw = y^^e-^erfc (-ayfijf) . 

as n — > oo, uniformly with respect to z in a neighborhood of z = 1. 

This expansion is in agreement with the general case 2-^1 (a, 6 — A; c + A; —2) considered 
in Olde Daalhuis (2002). 

The relation with the Jacobi polynomials is: 



2-f 1 



— n, 1 
n + 2 



. _,\. _ ( n + 1 ) ! fl , r p(n+l,-n-l) f IZLf 



See Figure 4 for the distribution of the zeros of 



P^iz), n = 30, a = n + l, /3 = -n-l. 
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• • • 



0.5 



-0.5 



Figure 4. Zeros of P 7 



(z), n = 30, a 



P = n+1. 



n 



As n — > oo, the zeros approach the curve |l — z 2 \ = 1. 
6. Asymptotics of some 3 F 2 polynomials 

In the previous examples we used integral representations of the Gauss function. The dif- 
ferential equation can also be used for obtaining asymptotic expansions. As a rule, methods 
based on differential equations provide more information on the coefficients and remainders 
in the expansions than methods based on integral representations. 

For the generalized hyper geometric functions asymptotic expansion are known when the 
argument z is large, and the other parameters are fixed; see Braaksma (1963) and Luke 



For the p F q — functions a differential equation is available of order max(p, q + 1). However, 
for higher order equations no methods are available for deriving (uniform) expansions for 
large values of parameters. 

The pF q — functions can be written as a Merlin-Barnes contour integral, but, again, no 
methods are available for deriving (uniform) expansions for large values of parameters from 
these integrals. 

The terminating p F q — functions (one of the parameters Ofc is equal to a non-positive in- 
teger) and of unit argument (z = 1) are of great interest in special function theory and in 
applications. 

In many cases recursion relations are available. For the examples to be considered here 
we use ad hoc methods. It is of interest to investigate recursion methods for the same and 
other problems. 

6.1. A first example 

As a first example (Larcombe & Koornwinder (2000)), consider the asymptotics of 



(1969). 
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It was conjectured (by Larcombe) that lim^-^ f(n) = 2 and Koornwinder gave a proof. By- 
using 

(*)fc = (~ 1 ) fcn! f 1 t k -hi-t) n - k -? dt 

(i-n) v^r(n +i )y (1 tj 
fc = 0, 1, . . . , n, we can write f(n) as an integral 

n! Z" 1 i, i ^ ( - n A t \ 

and the Gauss function is a Legendre polynomial. We obtain 

= m m / sin " 0Pn[^-Ade. (6.i) 

*\\)n(\) n Jo \sm0J 

There is no straightforward way to obtain asymptotics out of this. A few manipulations with 
(6.1) give the result 

n\n\n\ 

f(n) 



2"(i)n(i)n ^ 

where c n is the coefficient in the Bessel function expansion 

2 00 

e w ' 2 I (w/2)\ =Y J ^w n , 

n=0 

and now the asymptotics easily follows. We have 

1 f [e™/ 2 I (w/2)]\ 1 f e 2w \- W /2 T( 



and the contour C is a circle around the origin. The main contribution comes from the saddle 
point of J„ +1 , that is from w = wo = n/2. 
Because f(w ) ~ — , we have 



2 2" 

c n ~ r , n — »■ oo. 

7rn n! 

This gives finally f(n) ~ 2. A complete asymptotic expansion follows by using more terms 
in the expansion of f(w) at w = Wq. 

6.2. Generalizing Kummer's identity 

Kummer's identity reads: 



F ( a,b \ _ r(l + Q-b)r(l + iq) 

2 1 \,l + a-6' J r(l + a)r(l + ia-6) 



where 1 + a - b ^ 0, -1, -2, . . . 
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A generalization (Vidunas (2001)) is written in the form 



i + a-6' y w r(a)r(i a + i-6) ^ v ; r(a)r(i a -6)' 

Vidunas showed that, for n = —1, 0, 1, 2, ... , (a) n / and a — 6 / 0, — 1, —2, . . .: 
Q(n) = — 3 F 2 |^ ~ ^ + k ~h k + 5 " 6 • 1^) = - S F 2 ( + 5. " K & • i 



(6.2) 



We are interested in the asymptotic behaviour of P(n) and Q(n) for large values of n. 
Using, if $td > 5ftc> 1, 

^ = TV / (d ? TV ; f l 1?- c - l (l-t)<+ k - l dt, 
(d) k T(d-c)T{c)J K ' 

and a few manipulations of the Gauss functions, we obtain the integral representations, for 
> mb > 0, 

P{n) = ~T\ / shin 1 " 2 "" 1 (t) cosh"""^) cosh(n + l)tdt, 

r(6)r( 3 a- fe) y f6g . 

Q ^ = rmrfi x' M / sinha " 26_1 ^) cosh— (t) sinh(n+ l)t dt, 

We can use standard methods for obtaining the asymptotic behaviour of these integrals. First 
we write the hyperbolic functions cosh(n + l)t and sinh(n + l)t as exponential functions. 
We obtain for n — > oo 

r(6)T(2a-6)^ n b+ 2 fc 

2^T(|a+j) - r(b+ifc) 
r(6)T(2a+ 2 -6) ^ n b +2 fc 

where 

c = 1, ci = 0, c 2 = | (56 — 4a + 3). 

These expansions are not valid for 6 = ^a, 6 = \a + |, respectively. In these cases we 
have ^ 

as follows from (6.2). Similar for the cases 

b=-a + m, b=-a + - + m, m = 0,1,2,.... 
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For negative values of n different integral representations should be used. We have 



P( _ n _ i) = 2" a "t )ra 3^ ( ~^ ~M + ^J a " 6 ; l i . 



(1 " b), 

in (2 ~~ 2 a )™ 



Q(-n-l) =-n2" V l ^ 



(1 - b) n 

with the integral representations 

2 n+1 T(l - b) 



P(-n-l)= — f 
Q(-n-l) = 



r( i a _ 6)r( i_ i o) y 

2 n+l _ fe j 



r(i a + i-6)r(i-ia) J 



2 ' 2° n 

|n+ ±,-§n + l,±a + \ -b_ x 



cos n+1_a t sin 2c_1 t cos ntdt, 

n+1 " a t sin 2 "" 1 * sin ntdt. 



cos 



By integrating on contours (ioo,0)U(0, |7r)U(|7r, ^n+ioo), similar integrals arise as in (6.3), 
and again standard methods can be used for obtaining the asymptotic behaviour. 

7. Concluding remarks 

[1.] The asymptotic analysis of all 26 cases in 

/ a + eiA,fe + e 2 A \ 
2F1 ( „ 1 „ a ; z) , ej=0,±l, 



c + e 3 A 



can be reduced to the 4 cases 





ei 


e 2 


e3 


A 








+ 


B 







+ 


C 


+ 







D 


+ 


2+ 






[2.] These cases are of interest for orthogonal polynomials and special functions. 

[3.] New recent results on uniform asymptotic expansions have been published or an- 
nounced for special cases. 

[4.] The distribution of the zeros of Jacobi polynomials for non-classical values of the 
parameters a and f3 shows interesting features. New research is needed for describing 
the asymptotics of these distributions. 

[5.] The asymptotics of 3F2 terminating functions with large parameters is quite difficult. 
Standard methods based on integrals and differential equations are not available. Re- 
cursion relations may be explored further. 
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